Mesoscopic cross— film cryotrons: Vortex trapping and 
dc— Josephson— like oscillations of the critical current 
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We investigate theoretically and experimentally the transport properties of a plain Al supercon- 
ducting strip in the presence of a single straight current-carrying wire, oriented perpendicular to 
the superconducting strip. It is well known that the critical current of the superconducting strip, 
I c , in such cryotron-like system can be tuned by changing the current in the control wire, I w . We 
demonstrated that the discrete change in the number of the pinned vortices/antivortices inside the 
narrow and long strip nearby the current-carrying wire results in a peculiar oscillatory dependence 
of I c on I w . 

PACS numbers: 74.25. F- 74.25.Sv 74.78.-w 74.78.Na 



I. INTRODUCTION 

The original idea to control the resistance of a long 
superconducting (S) wire by means of a magnetic field 
generated by a coil wound locally over the wire, was 
proposed by Buck in 1956. 1 If the magnetic field in- 
side the solenoid with permanent driving current exceeds 
the critical field of the central type-I superconducting 
wire, superconductivity will be completely suppressed 
and this device (cryotron) will be switched from a low- 
resistive state to a high-resistive state. Further investiga- 
tions in the 1960'a 2 - - — revealed that the cryotrons can be 
potentially used as superconducting computer elements 
(switches, binary adders, shift registers, AND/OR gates 
etc). However the performance of such elements at high 
frequencies was found to be worse than ordinary semicon- 
ducting chips, since the cryotrons were rather slow and 
energy-consuming. Furthermore, the miniaturization of 
cryotron-like devices, which appears to be important for 
reducing the characteristic time constant, was very lim- 
ited at those times. 

Revival of interest to superconductor-electromagnet 
(S/Em) hybrids came in the 1990's in connection with 
the problem of the influence of a spatially-modulated 
magnetic field on superconductivity. The development 
of advanced techniques for material deposition and litho- 
graphic methods have made it possible to fabricate com- 
posite structures with controlled arrangement of super- 
conducting, normal metallic, ferromagnetic and insulat- 
ing layers at micron and submicron levels. Panncticr 
et al& demonstrated that the dependence of the critical 
temperature T c on the applied magnetic field H for such 
a hybrid system consisting of a plain Al film and a litho- 
graphically defined array of parallel metallic lines can be 
nonlinear and even non-monotonous in contrast to a plain 
superconducting film in a uniform magnetic fields Such 
tunability of the T C (H) dependence by the stray field of 



the meander-like wire reflects directly the controllable 
modification of the standard Landau spectrum for the 
order parameter (OP) wave function in a periodic mag- 
netic field£~— It is worth noting that keeping the nonuni- 
form magnetic field of the control coils/wires requires ex- 
penditure of energy and may result in parasitic heating 
effects. Therefore, later on the interest was turned to 
superconductor-ferromagnet (S/F) hybrids since the in- 
homogeneous magnetic field in S/F systems, conditioned 
by the nonuniform distribution of magnetization, can be 
obtained without energy costs. The influence of a nonuni- 
form magnetic field produced by ferromagnetic elements 
on nucleation of superconductivity and low-temperature 
properties of superconducting specimens was studied in- 
tensively for flux-coupled S/F hybrids during the last two 



decades (see review: 



,11-13 



and references therein). We 



would like to note that the amplitude of the stray mag- 
netic field and its profile are dictated by the saturated 
magnetization of the fcrromagnet and by the shape of 
the ferromagnetic elements, which eventually restricts 
the flexibility of the S/F hybrids. 

In the present paper we report on a so far unexplored 
limit of a cryotron-like system at micro- and nanoscales 
bringing together the ideas and approaches developed for 
both S/Em and S/F hybrids. With that purpose we fab- 
ricated a composite structure consisting of a type-II su- 
perconducting strip and a single straight current-carrying 
S wire oriented perpendicular to the strip. Since the stray 
magnetic field of the straight wire is maximal near the 
wire and it decays approximately as 1 jr at large distances 
r from the wire, one can expect that only a small part of 
the superconducting strip in close vicinity to the wire will 
be affected by the nonuniform magnetic field.— A simple 
change in the current in the wire, I w , gives us the op- 
portunity to control the maximal value of the magnetic 
field, Bo, generated by the current in the control wire. 
By varying Bq one can create different states with par- 
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tially and completely suppressed superconductivity in the 
considered cross-film cryotron. Thus, the main point of 
our research is to "scan" all possible intermediate states 
lying between uniform superconductivity (for small -Bo) 
and fully depleted superconductivity (for large Bq) and to 
study the effect of the transitions between various super- 
conducting states on the critical current for mesoscopic 
type-II cross-film cryotrons. 



II. THEORY 
A. Model 

To simulate the onset of vortex dynamics in the super- 
conducting strip in the field of the current-carrying wire 
(i.e. in the cross-film cryotron) under the action of the 
transport current we use a time-dependent Ginzburg- 
Landau model.— The thickness D s of the superconduct- 
ing strip is assumed to be infinitely small, so the effec- 
tive magnetic field penetration length A = \\/D s ex- 
ceeds substantially the lateral dimensions of the super- 
conducting sample (Al is the London penetration depth). 
Consequently, both the magnetic field B = rot A and 
the vector potential A are determined solely by external 
sources. We use the following units: m*<r n (3/(2e 2 a) for 




FIG. 1: (color online) (a) Schematic presentation of the meso- 
scopic cross-film cryotron, (b) Profile of the z— component of 
the magnetic field, induced by the current-carrying wire, in 
the middle plane of the superconducting strip. The maximal 
and minimal b~ values correspond to the edges of the current- 
carrying wire (shown in grey). 



time, the coherence length £(0) at temperature T = 
for distances, $o/(27r^(0)) for the vector potential and 
4eci 2 £(0)/ (Hj3) for the current density, where a = -or 
and /? are the conventional parameters of the Ginzburg- 
Landau expansion, r = (1 — T/T c q), T c $ is the critical 
temperature at B = 0, e and m* are charge and the 
effective mass of carriers, er„ is the normal state con- 
ductivity. Then the time-dependent Ginzburg-Landau 
equations take the formic 

u (jt + Up} i/> = r - |V>| V) + (V + zA) 2 V, (f ) 

VV = divj s , j s = -^t{v>*(V + zA)?A-c.c.}. (2) 

Here ijj is the normalized order parameter (OP), tp is the 
dimcnsionless electrical potential, j s is the dimensionless 
density of superconducting currents, u is the rate of the 
OP relaxation, c.c. stands for complex conjugate. To 
complete the problem we apply the boundary conditions 

(s + "-) r *-°- (l) r = 3 "" (3) 

which seems to be natural for the superconduc- 
tor/vacuum or superconductor/insulator interfaces, n is 
the normal vector to the sample's boundary T, j ext is the 
normal component of the inward (outward) flow of the 
external current. To detect a vortex on the grid cell cor- 
responding to the certain grid node (xi,yj), we use the 
following criterion^ 

+ arg(VCf u+iV'i+ij) + axg(ipt + ijipi,j) = ^N, (4) 

i-'i.j = i>(xi,yj), ar g(V') is the argument of the complex 
number ipi j, belonging to the interval (-7T, it]. If N = +1, 
there is a vortex in the cell, while N = — 1 corresponds 
to an antivortcx. It is important to note that we define 
the vortex (antivortex) as a position of the singularity 
of the OP phase, 9 = arg(?/>), however the direction, in 
which the superconducting currents circulate around the 
center of the vortex (antivortex), depends on the charge 
of carriers (i.e. the e sign). 

To facilitate further comparison between theory and 
experiment (Section III), for our modeling we choose 
the parameters, typical for mesoscopic Al-based super- 
conductors: the coherence length £(0) = 0.15 fj,m at 
T = 0, length and width of the strip L x = 9 /im and 
L y = 3 /im, width and thickness of the current-carrying 
wire L w = 1.5 /zm and D w = 0.05 /im, and the separation 
between the strip and the wire h = 0.05 /im (Fig. |T|) . 

B. Mixed state in cross-film cryotron in 
equilibrium at H = 

The stray field generated by the control wire is uniform 
across the superconducting bridge (in the y— direction, 
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FIG. 2: (color online) Contour plots of \ip\, showing equi- 
librium vortex-antivortex patterns in the thin superconduct- 
ing strip (60 ^(O) x20 ^(O) in size), appearing at H = and 
T — 0.9 T C Q in the field of the current-carrying wire: I w = 5 
mA (a), I w = 8 mA (b), 7 TO = 10 mA (c), and 7 TO = 14 mA (d). 
Lighter shades correspond to higher \ty\ values, darker shades 
are the regions with suppressed superconductivity. The grey 
rectangle in the center of the strip is the projection of the 
current-carrying wire. The symbols depict the OP phase sin- 
gularities (i.e., the cores of the vortices and antivortices). 



Fig.rjJ)). However in the a;— direction the field is substan- 
tially inhomogeneous with the 6 z -component reaching its 
extreme values near the edges of the current-carrying 
wire. The increase in the control current I w leads to 
a gradual enhancement of the maximal b z value. This 
eventually results in the local suppression of the energy 
barrier for vortex penetration from the edges of the strip 
near the \b z \ maxima. These regions play the role of pre- 
defined gates for vortex entry. 

The stray field of the control wire acts as (i) a tun- 
able source of vortices and antivortices and, (ii) a vortex 
trap which prevents vortex-antivortex pairs from their 
mutual annihilation and the escape from the supercon- 
ducting bridge. The latter follows from the distribution 
of the screening superconducting currents j s induced by 
the current-carrying wire (see Fig. \5jp). Indeed, the dif- 
ferent signs of the x— component of j s near the top and 



bottom edges of the strip mean that the Lorentz force 
f*L = ($o/c) [js x z ] acting on the vortex (antivortex) 
near these edges will be oriented in such a way in or- 
der to push the vortex (antivortex) to the inner part of 
the strip. In the absence of an external magnetic field 
H the appearance of a vortex should be accompanied by 
the formation of a symmetrically positioned antivortex 
and the number of vortex-antivortex pairs N increases 
as \I W \ increases (Fig. [2]). Since the b z — field roughly 
decays inversely proportional to the distance from the 
wire, the vortices and antivortices are also distributed 
non-uniformly: the closer to the wire, the smaller the 
distance between two vortices (antivortices) is. 
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FIG. 3: (color online) (a) Critical current I c as a function of 
control current I w calculated at H = and T/T c o = 0, 0.5, 
and 0.9. For each segment of the I c — I w diagram we indi- 
cate the number of the trapped vortices/ antivortices JV in the 
strip, (b-c) Left panels show the stationary OP wave func- 
tions \ip\ calculated for T/T c n = 0.9 and I m = 3.4 mA for 
two different values of the injected current [points A and B in 
panel (a)] . The grey rectangle in the center of the strip is the 
projection of the current-carrying wire. The correspondent 
right panels show the profiles of the order parameter wave 
function \ip\, the longitudinal x— components of the current 
j 3tX and the supercurrent velocity v s<x oc (80/dx — 2ttA x /<&q) 
across the strip (i.e. along the dashed line). 
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C. Critical current of cross-film cryotron at H — 

Depending on the magnitude of the bias current I ex ti 
injected into the superconducting strip, there are two 
distinct regimes. After applying a small bias current 
the disturbed vortex and antivortex tend to relax into a 
new stationary configuration. This state with motionless 
vortcx-antivortex pairs is characterized by the absence of 
an electrical field inside the superconductor (except for 
tiny regions near the boundaries where the current in- 
jection takes place) and by a zero voltage drop. If I ext 
exceeds the critical current I c , the stability of the vortex- 
antivortex ensemble breaks down: vortex and antivortex 
periodically in time enter and exit into/from the super- 
conductor. This means that the sample has switched to 
a dissipative regime. The calculated dependencies I C (I W ) 
at H = for different temperatures arc shown in Fig. [3^. 
We see that I c decays with oscillations as I w sweeps and 
the I c oscillations become more pronounced at high tem- 
peratures. 
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FIG. 4: (a) The calculated dependencies I C (I W ) for T/T c o = 
0.9 and H = (black line) and H/H c2 = 0.1 (blue line). For 
each segment of the I c — Iw diagram we indicate the number 
of vortices N v and antivortices N av , trapped in the supercon- 
ducting strip, (b-e) Examples of the stationary OP patterns 
at T/T c o = 0.9 and H/H C 2 = 0.1 and at bias current close to 
the corresponding critical values. 

We can explain the appearance of the I C {I W ) oscilla- 
tions as follows. If H = and I w = 0, superconductivity 
will survive until the injected current density exceeds the 
depairing limit. For nonzero I w the transition from the 



non-dissipative state to the resistive regime occurs via 
permanent formation of vortex-antivortex pairs at the 
opposite edges of the superconducting bridge and their 
motion across the bridge (see inset in Fig. [3]). Due to the 
superposition of the injected current and the currents, 
induced by the wire, the limit for the vortex-antivortex 
pair generation can be reached at smaller values of the 
bias current; therefore I c monotonously decreases as /„, 
increases. A further increase in I w leads to (i) a decrease 
in the threshold value of the bias current required for the 
creation of the first vortex, and (ii) an enhancement of 
the trapping potential for vortices. Both consequences 
make it possible to stabilize the vortex-antivortex pair in 
the presence of the transport current (Fig. (3^). The self- 
currents generated by the vortex and antivortex partly 
compensate the superflow conditioned by the bias cur- 
rent near the regions with suppressed OP wave function. 
Therefore these "gates" will be closed and the formation 
of new vortices and antivortices will be impeded unless 
the maximal value for supercurrent velocity v s , which 
is proportional to V(9 — 27rA/<I>o, again reaches its crit- 
ical value at the strip edges. This means that after the 
first vortex-antivortex pair is trapped, the critical current 
will be larger than that in the Meissncr state. Then this 
process is repeated periodically as I w increases, result- 
ing in periodic variations in I C (I W ). Thus, the discrete 
change in the number of the pinned vortices and antivor- 
tices in the superconducting strip of a finite width can be 
identified by considering the position of the cusps on the 
Ic(Iw) curve. In contrast to the I c oscillations observed 
in mesoscopic superconducting squares in a uniform mag- 
netic fields the nonuniform magnetic field gives us the 
possibility to see the quantization effects in a long super- 
conducting strip, due to the confinement potential along 
the strip produced by the control wire. 

The absolute value of the magnetic flux $!/ 2 piercing 
a half of the superconducting strip, 



Lx/2 L y /2 



$1/2 



\b z (x,y)\ dxdy, 



(5) 



-Ly/2 



is a linear function of I w . In a certain sense the oscil- 
lations of the critical current of the cryotron as I w (or 
$1/2) varies, are similar to the standard Fraunhofer pat- 
tern describing the dependence of the critical current of a 
Joscphson junction on the total magnetic flux <f> through 
the junction area. 19 ' 20 The built-in field of the wire guar- 
antees reversible entrance of vortices and antivortices into 
the superconducting strip as I w increases, and the tun- 
able modification of the distribution of the OP phase in 
the restricted part of the strip by the trapped vortices 
and antivortices. This area with partly suppressed super- 
conductivity near the control wire acts like an effective 
weak link, since the vortex dynamics in this area deter- 
mines the flow of the bias current and the resistance of 
the entire strip. 
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D. Mixed state and critical current of cross-film 
cryotron at H ^ 

The external magnetic field H, applied perpendicularly 
to the sample's plane, breaks the symmetry between vor- 
tex and antivortex since the total magnetic flux piercing 
the sample is not zero anymore. As a result, the entrance 
of a new vortex is not generally accompanied by the en- 
trance of an antivortex: the numbers of trapped vortices 
and antivortices, N v and N av , may differ in contrast to 
the previously considered case H = 0. 

The results of the calculations for T/T c q = 0.9 and at 
H/H C 2 = and H/ H c i = 0.1 are compared in Fig.|4ji. It 
is clear that the oscillating behavior in the I c vs. I w re- 
mains even in the presence of the applied magnetic field. 
Since the constructive superposition of (i) the injected 
current, (ii) the currents induced by the wire and (iii) 
the currents induced by the external field, enhances the 
local flow of superconducting condensate, the conditions 
for vortex entry will be fulfilled at smaller I w values. This 
results in a shift of all the cusps on the I C {I W ) dependence 
towards to zero. The non-equivalence between vortices 
and antivortices allows us to get new stable states of re- 
duced symmetry (Fig. 2]3-e) which are forbidden for the 
symmetrical sample at H = 0. The transition between 
such states sweeping I w are accompanied by additional 
cusps, which leads to more complex oscillatory properties 
of the I c on I w dependence. 



III. EXPERIMENT AND DISCUSSION 

To study the transport properties of the mesoscopic 
cross-film cryotron, we fabricated the hybrid structures 
consisting of a 4 ^m wide and 120 nm thick super- 
conducting Al strip on top of a 1.5 fxm wide and 50 
nm thick current-carrying Nb wire. The Nb thin film 
was deposited via dc magnetron sputtering at a rate of 
2.5 nm/s on an oxidized Si substrate at room tempera- 
ture under the pressure of Ar plasma of 6xl0 -3 mbar. 
The Nb wire was fabricated by e-beam lithography in 
combination with argon ion milling. The Al film was 
grown in a molecular-beam epitaxy apparatus at a rate 
of 0.2 nm/s at room temperature and nominal pressure of 
1 x 10 -8 mbar. The Al strip was patterned by standard 
e-beam lithography, subsequent evaporation and lift-off 
technique. A 120 nm thick Ge layer separates the two 
metals from each other to avoid electrical contact be- 
tween them. We tested more than 30 identical hybrid 
samples prepared at different times but under similar 
conditions. At least 10 samples have no problems with 
electrical leakage, since the applying a voltage between Al 
strip and Nb wire leads to an absence of detectable cur- 
rent. It gives us a lower limit on the contact resistance of 
the order of 5 Mf2. An Atomic Force Microscopy (AFM) 
image of the sample is shown in Fig. [5j 

Considering the dependence of the sample resistance R 
on H and its displacement upon varying T , one can com- 




FIG. 5: (color online) AFM image of the sample: the top 
(light) element is the superconducting Al strip with two per- 
pendicular contacts for the measurement of the voltage drop, 
the bottom element oriented perpendicular to the strip is the 
current-carrying Nb wire. 



pose the phase transition line T C (H) shown in Fig. ^p. 
By identifying this line with the temperature depen- 
dence of the critical field of surface superconductivity 
Has = 1.69 (1-T/T c0 ), we estimated T c0 ~ 1.265 K, 
the upper critical field ~ 104 Oe and the coherence 

length £(0) = (^q/2ttH^° ) ) 1/2 ~ 175 nm extrapolated to 
T = 0. The mean free path I ~ 35 nm can be obtained 
from the value of the normal resistivity p n at low temper- 
atures according to the formula:— p n I ~ 6 x 10~ 16 fi-m 2 . 
The i ~ 35 nm value gives us the London penetration 
depth and the effective magnetic field penetration depth 
A L (0) ~ 175 nm and A(0) = A|(0)/d ~ 250 nm re- 
spectively in dirty limit at T — 0. All obtained values 
are in accordance with previous studies on thin film Al 
structures^ 

Figure [7K shows the I c of the Al strip as a function of 
the current I w in the control Nb wire, measured at H = 
and different temperatures close to T c q. Experimentally 
the critical current was determined from the isothermal 
current (J) - voltage (V) dependencies for different volt- 
age criteria (from 0.1 fiV up to 10 pV between the two 
voltage contacts at a distance 20 micron) for the mea- 
sured voltage drop.— For our measurements the accu- 
racy SI C of the determination of the critical current was 
3 yLtA. This means that within the interval 2 5 I c we ob- 
served a sharp transition from the superconducting state 
(typical noise level of the order of 5 x 10 -8 V) to the 
resistive state. The plateau in I C (I W ) might be associ- 
ated with the influence of sample imperfections, which 
facilitate the formation of vortices near defects even in 
zero magnetic field at current densities smaller than the 
theoretical limit. In accordance with our simulations 
the monotonously decreasing part of I C (I W ) has to be 



attributed to the transition from the Meissner state to 
the resistive state with a single moving vortex-antivortex 
pair. The next stages correspond to the stabilization of 
the first and all subsequent vortex-antivortex pairs in the 
sample. The switching between these states at varying 
I w corresponds to the local I c minima. We clearly ob- 
served the reproducible oscillations of the critical current 
for all measured samples. 

The effect of the external field is illustrated by Fig.[7j3. 
We observed the anticipated shift of the monotonously 
decreasing part in I C (I W ) towards smaller I w values, con- 
ditioned by the constructive superposition of the injected 
current, the currents induced by the wire and by the ex- 
ternal field. This ensures that the conditions for the first 
vortex entry will be fulfilled at smaller I w values. It is in- 
teresting to note that the shape of the oscillating I C (I W ) 
curves becomes more complicated as H increases. We 
suppose that the external-ficld-induccd modification of 
I C (I W ) reflects the formation of exotic non-symmetrical 
vortex states and their depinning under the action of the 
bias current. 

To compare theory and experiment, we plot the posi- 
tions of the minima and maxima in the I C (I W ) depen- 
dence at H — (Fig. [5J. It is easy to see that our 
simple model describes both the general oscillating be- 
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FIG. 6: (color online) (a) Typical dependence of the voltage 
drop, V, as a function of H at I w = for different temper- 
atures. The bias dc current I is equal to 50 ^A, the normal 
state resistance R n — 0.45 Q. (b) The phase boundary T C (H) 
extracted from the magnetoresistive measurements at I w — 
according to the criterion R(H,T) = 0.99 R n . The dashed 



line is the linear approximation (1 — T/T c q) = 0.59 \H\/H : 
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104 Oe, T c0 = 1.265 K. 



FIG. 7: (color online) (a) Measured dependence 7 C on I w at 
H = for different temperatures, (c) Measured dependence 
I c on I w at T — 1.16 K for different H values. Both insets 
in the panels (b) and (c) show zoomed parts of the I C (I W ) 
curves (within the box) in semi-logarithmic scale. The arrows 
in the inset (a) show the positions of the extremal 7 C values 
for T = 1.16 K and H = 0, which are compared with theory 
in Fig. [8] The confidence intervals are smaller than the size 
of symbols, and therefore error bars cannot be shown in this 
scale. 



havior for the I C {I W ) dependence for real sample and the 
observed period AI W quite well. The fact that all theo- 
retical points in Fig. [8] lie below the experimental data, 
may be caused, e.g., by a nonuniform current distribu- 
tion in the control wire, or by a small difference between 
the parameters of the model problem and the dimensions 
of the tested samples. Nevertheless using the estimated 
period AI W ~ 1.45 mA, one can calculate the change in 
the magnetic flux A^ 1 / 2 through the half of the super- 
conducting strip and corresponding to the appearance of 
the vortex-antivortex pair: A$i/2 — l-l'&o- 

In this paper we have reported on the first results ob- 
tained for the prepared S/Em samples. We note that 
the oscillatory behavior of the critical current as a func- 
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FIG. 8: (color online) Positions of the I c minima (left panel) 
and the I c maxima (right panel), calculated at H = and 
T/T c q = 0.9 (red circles) and measured at T — 1.16 K 
(T/T c o ~ 0.91, black dots). Dashed red lines are guide to 
the eyes. The positions of the maxima and minima deter- 
mined experimentally are also marked by vertical arrows in 
the inset of Fig. [7^. 



carrying wire, oriented perpendicular to the strip. The 
stray field, generated by the control wire, makes it pos- 
sible to tune the transport properties of the long super- 
conducting strip locally by means of creation and pin- 
ning of vortices and antivortices. This area, where vor- 
tices and antivortices become confined, plays the role of 
a "bottleneck" for the transport current, since the mo- 
tion of vortices and antivortices is the main source of the 
dissipation. Thus, the peculiar configurations of vortex- 
antivortex pairs determine the critical current I c of the 
S strip. The transition between different stable vortex 
patterns upon varying I w results in an oscillatory de- 
pendence of I c on I w . The observed effect, which looks 
similar to the Fraunhofer oscillating pattern for conven- 
tional Joscphson junctions, can be potentially interesting 
for the design of superconducting interferometer devices 
based on S/F or S/Em hybrid structures. 



tion of I w was found to be reproducible for all measured 
samples, which had identical geometry and size as the 
presented device. However the issues concerning the in- 
fluence of the dimensions of the superconducting strip 
and the current-carrying wire on the I c oscillations re- 
main unaddressed. 

IV. CONCLUSION 

We studied the transport properties of a hybrid sys- 
tem consisting of a superconducting strip and a currcnt- 
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